Electron scattering in isotonic chains as a probe of the proton shell structure of 

unstable nuclei 



X. Roca-Maza 1 '^ M. Centelles 1 , F. Salvat 1 , and X. Vinas 1 
1 Departament d'Estructura i Constituents de la Materia and Institut de Ciencies del Cosmos, 
Facultat de Fisica, Universitat de Barcelona, Diagonal 645, 08028 Barcelona, Spain 
2 INFN, sezione di Milano, via Celoria 16, 1-20133 Milano, Italy 

Electron scattering on unstable nuclei is planned in future facilities of the GSI and RIKEN up- 
grades. Motivated by this fact, we study theoretical predictions for elastic electron scattering in 
the N = 82, N = 50, and TV = 14 isotonic chains from very proton-deficient to very proton-rich 
isotones. We compute the scattering observables by performing Dirac partial-wave calculations. The 
charge density of the nucleus is obtained with a covariant nuclear mean-field model that accounts 
for the low-energy electromagnetic structure of the nucleon. For the discussion of the dependence 
of scattering observables at low-momentum transfer on the gross properties of the charge density, 
we fit Helm model distributions to the self-consistent mean-field densities. We find that the changes 
shown by the electric charge form factor along each isotonic chain are strongly correlated with the 
underlying proton shell structure of the isotones. We conclude that elastic electron scattering ex- 
periments in isotones can provide valuable information about the filling order and occupation of the 
single-particle levels of protons. 
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I. INTRODUCTION 

Since the 1960's, elastic electron scattering has been 
utilized to obtain accurate information on the size and 
shape of nuclei . Because electrons and nucleons in- 
teract essentially through the electromagnetic force, the 
nucleus remains rather unperturbed during the scatter- 
ing process and the analysis of the data is not hampered 
by uncertainties associated with the strong interaction. 
Thus, electron scattering is able to provide very clean 
information about the charge distribution of atomic nu- 
clei 

Low-energy nuclear physics is nowadays moving very 
fast towards the domain of exotic nuclei [9[ . This is due to 
the development of successive generations of radioactive- 
isotope beam (RIB) facilities p^ - fTEj . such as FAIR 
and SPIRAL2 in Europe, FRIB in North America, and 
HIRFL-CSR, RARF or RIBF in Asia, which will allow 
studying the properties of nuclei beyond the stability val- 
ley. Many interesting effects have already been discov- 
ered in exotic nuclei, such as neutron and proton halos, 
neutron skins, and new magic numbers. These effects 
may be related to the structure of the nucleon distri- 
butions far from stability. As with stable nuclei, one 
way of exploring the structure of exotic nuclei is through 
the electromagnetic interaction. For this purpose, a new 
generation of electron- RIB colliders using storage rings 
is under construction by RIKEN (Japan) [lj| [l6| and 
at GSI (Germany) [H |18| . It is expected that in the 
near future the SCRIT project in Japan[l9|-[2l| and the 
ELISe experiment at FAIR in Germany (22[ will offer the 
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opportunity of studying the structure of unstable exotic 
nuclei by means of electron scattering. 

On the theoretical side, much work has been devoted to 
the study of charge distributions of exotic nuclei through 
calculations of both electron scattering (see e.g. Refs. I23t- 
I29I ]) and proton scattering (see e.g. Ref. [3(|). Suda |3l| 
pointed out that in electron scattering off unstable nuclei 
the maxima and the minima of the charge form factor are 
very sensitive to the size and the diffuseness of the charge 
density. This fact has been confirmed by different works 
that have analyzed the behavior of the char ge f orm factor 
along isotopic [23|, H3, H3, [H| and isotonic [29| ] chains. 

To probe the charge distribution in nuclei, the electron 
beam energy needs to be of the order of a few hundred 
MeV. As one deals with relativistic electrons, it is manda- 
tory to solve the elastic scattering problem of Dirac par- 
ticles in the potential generated by the nuclear charge 
density. The simplest approach is the plane-wave Born 
approximation (PWBA) where the initial and final states 
of the electron are described by Dirac plane waves. The 
PWBA accounts for many features of electron scattering 
but it cannot provide an accurate description of the elec- 
tric charge form factor, in particular near the deeps. The 
most elaborated calculations of electron-nucleus scatter- 
ing are obtained by the exact phase-shift analysis of the 
Dirac equation. This calculation scheme is known as 
distorted-wave Born approximation (DWBA) [32J and 
has been used to analyze different aspects of the scat- 
tering of electrons by nuclei, see e.g. Refs. I24I l28l. l33l| 
and references therein. It is worth noting that for high- 
energy electron scattering, the eikonal approximation is 
a reliable choice to deal with the problem (34J; it has 
been applied in [27], H^, [35[ to study electron scattering 
off proton- and neutron-rich nuclei. 

The charge density of the target nucleus is one of the 
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basic ingredients of the electron-nucleus scattering prob- 
lem. For medium and heavy nuclei, the theoretical charge 
densities can be calculated in the mean-field approxima- 
tion using non-relativistic nuclear forces or relativistic 
mean field (RMF) models. It is known that the over- 
all trends of the elastic scattering of electrons by stable 
medium and heavy nuclei, are well reproduced by the 
mean-field charge densities computed with nuclear mod- 
els that have been calibrated to describe the ground-state 
properties (in particular the charge radii) of some se- 
lected nuclei. However, different nuclear models differ in 
the fine details and describe with different quality the 
experimental scattering data. See Ref. [28| for a recent 
comparison of the elastic electron scattering results pre- 
dicted by different nuclear mean-field models. 

In Ref. [28[ we studied elastic electron scattering along 
the Ca and Sn isotopic chains in DWBA. In that work 
we reported several correlations among scattering observ- 
ables and some properties of the nuclear charge density 
along the isotopic chains [28|. In the present work, we 
investigate what information on nuclear structure can be 
gained from the study of elastic electron scattering in 
the N — 82, 50, and 14 isotonic chains, which cover dif- 
ferent regions of the mass table. We aim at extracting 
general trends, according to current mean-field theories, 
about the behavior of some observables that may be mea- 
sured in experiments performed with unstable nuclei in 
the low-momentum transfer region. Our choice of these 
isotonic chains is motivated by the following reasons. On 
the one hand, information about the structure of unsta- 
ble nuclei belonging to the N — 82 and N — 50 isotonic 
chains may be relevant in the astrophysical context be- 
cause some of these nuclei could correspond to waiting 
points in r-processes [36|, [37} • On the other hand, scat- 
tering data for light nuclei, such as e.g. those of N = 14, 
are likely to be obtained in future electron scattering fa- 
cilities such as SCRIT (HHI and ELISe [Hj]. 

The study of elastic electron scattering along isotopic 
and isotonic chains explores different aspects of the nu- 
clear charge density. The electric charge form factor 
along an isotopic chain gives information about the effect 
of the different number of neutrons on the charge den- 
sity, which becomes more and more dilute and extends 
to larger distances as the neutron number increases [28| . 
In an isotonic chain, the changes in the charge form fac- 
tor primarily inform about the effect of the outer proton 
single-particle orbitals that are being filled as the atomic 
number increases in the chain. Thus, our previous [28j 
and present study together provide a survey of the evo- 
lution of the charge form factor with the neutron and 
proton numbers in different mass regions of the nuclear 
chart. 

The rest of this article is organized as follows. In Sec- 
tion II, we summarize the method employed in our study 
of electron scattering in isotonic chains. As the basic 
methodology follows that of Ref. [Hj|, we address the 
reader to that work and references therein for more de- 
tails about the relativistic nuclear mean-field theory and 



about the Dirac partial-wave analysis, which we perform 
using the ELSEPA code [38[ adapted to the nuclear prob- 
lem. We devote Section III to the presentation and anal- 
ysis of our numerical results for elastic electron-nucleus 
scattering in the TV = 82, 50, and 14 chains. Finally, our 
conclusions are laid in Section IV. 



II. METHOD 

To investigate electron scattering in isotonic chains we 
follow the method developed in Ref. 28]. For complete- 
ness, we summarize here the main aspects of this method. 
The electron beam energy in our investigation is fixed at 
500 MeV, which is a typical energy in electron-nucleus 
scattering experiments. Indeed, rather than discussing 
directly the differential cross section (DCS), we study 
the DWBA electric charge form factor because it is al- 
most independent of the electron beam energy in the low- 
momentum transfer regime, as it can be seen from Fig. 5 
of Ref. [28| and from the lower panel of Fig. @] below. We 
compute the electric charge form factor as follows: 

™i'-(b)(tt)"- w 

where da/dQ, is the DCS calculated in DWBA and 
da Mo tt/dfl is the Mott DCS. We use the exact Mott DCS 
(i.e., the DWBA for a point nucleus) instead of the DCS 
calculated in PWBA. Therefore, the charge form factor 
defined by Eq. ([1]) goes beyond the first Born approxi- 
mation. It will be denoted by -Fdwba(<7) hereinafter. 

We calculate the charge densities with the relativistic 
mean- field parametrization G2 (39ll40j |. which we also em- 
ployed in Ref. [Hj]. This nuclear model was constructed 
as an effective hadronic Lagrangian consistent with the 
symmetries of quantum chromodynamics. The nucleon 
density distributions are obtained self-consistently at the 
mean-field level by numerical solution of the correspond- 
ing variational equations. In contrast to most of the 
nuclear mean-field models that assume point densities, 
the G2 Lagrangian incorporates the low-energy electro- 
magnetic structure of the nucleon through vector-meson 
dominance [H, H(| ■ This implies that the charge density 
is obtained directly from the self-consistent solution of 
the mean-field equations without any extra folding with 
external single-nucleon form factors. It has been shown 
[39M4l| that the G2 relativistic mean-field interaction is 
a reliable parameter set both for calculations of ground- 
state properties of nuclei and for predictions of the nu- 
clear equation of state up to supra- normal densities, as 
well as for predictions of some properties of neutron stars. 
First calculations of the char ge f orm factor in PWBA 
with G2 were reported in Ref. |39| . 

In our calculations, we assume spherical symmetry. 
The suitability of this assumption for the N = 82, 50, 
and 14 isotones is confirmed, excluding a very few ex- 
ceptions of small deformation, by mass tables computed 
either with the non-relativistic Skyrme model [42j or with 
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the relativistic model [43|. Pairing correlations are im- 
portant for describing open-shell nuclei. We take pairing 
into account through a modified BCS approach that sim- 
ulates the continuum (needed for nuclei at the drip lines) 
through quasi-bound levels which are retained by their 
centrifugal barrier (neutron levels) or by their centrifugal- 
plus-Coulomb barrier (proton levels) The pairing 

interaction in this approach is described by means of a 
constant matrix element fitted to reproduce the exper- 
imental binding energies of some selected isotopic and 
isotonic chains as described in Ref. [44 |. It is to be men- 
tioned that a mean-field treatment is not expected to be 
sufficient for light exotic nuclei [25j. Thus, the N = 14 
isotonic chain studied below (with nuclei from 22 O to 
34 Ca) corresponds to a somewhat limiting case, and the 
mean-field results should be taken as semiquantitative. 
The calculations with the G2 model predict a relatively 
magic character of the neutron numbers N — 14 and 
N = 16. These neutron numbers have attracted some 
attention in recent theoretical and experimental studies 
as possible new magic numbers in exotic nuclei pBI - l49l |. 

The use of modeled charge densities and electric charge 
form factors in the experimental analysis of scattering 
data has been extensive in the past, and continues to 
date. This is because in many cases the parameters of 
the modeled charge densities are directly related with the 
size of the bulk and surface regions of the nucleus under 
study. In this way, the parametrized densities help to pro- 
vide a clear physical interpretation of the electron scat- 
tering data. This is the case of the so-called Helm model 
[50l | that we used for some calculations in our previous 
study of isotopic chains [28] . As one of our aims here is to 
predict global trends of the electric charge form factor in 
the low-momentum transfer regime in isotonic chains, we 
also will perform calculations with Helm model densities 
fitted to the self-consistent mean-field charge densities. 
We briefly summarize this procedure in the next subsec- 
tion. 



A. Equivalent Helm charge densities 

The original version of the Helm model j50j has been 
extended in various ways for a more accurate descrip- 
tion of the experimental charge densities (5lT[53j . In the 
simplest version of the model [5(J, the charge density is 
obtained from the convolution of a constant density po 
in a hard sphere of radius Rq with a Gaussian distribu- 
tion having variance a 2 . By construction, Rq gives the 
effective location of the position of the nuclear surface, 
and hence characterizes the size of the density profile, 
whereas the parameter a is a measure of the thickness of 
the surface region of the density distribution. The Helm 
charge density is then given by 



where 
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The two parameters, Rq and <r 2 , of the Helm model de- 
termine the charge density as well as the electric charge 
form factor within the PWBA: 
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where ji(x) is the spherical Bessel function. Note that 
we use natural units throughout the present paper. 

We proceed as suggested originally in Ref. [5fJ] to ob- 
tain the Helm parameters associated to a given nucleus. 
First, we require that the first zero of Eq. (j4j coincides 
with the first zero of the mean-field PWBA charge form 
factor (Fourier transform of the charge density obtained 
with the G2 model). We will refer to this charge form 
factor as Fpwba(</) hereinafter. Therefore, the radius of 
the equivalent Helm density reads 



Ro = — . 

qa 



(5) 



where x = 4.49341 is the first zero of ji(x) and qo is 
the momentum transfer corresponding to the first zero 
of Fpwba(<7)- Second, we determine the variance a 2 
of the Gaussian distribution such that |F( ff '(g max )| = 
|FpwBA(9max)|, where q max is the momentum transfer 
corresponding to the second maximum of |Fpwba(<z) 
(the first maximum appears always at q ~ fm). Us- 
ing Eq. one easily obtains 
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The PWBA form factor obtained in the Helm model 
neglects the distortion of the electron wave functions due 
to the electrostatic potential of the nuclear charge dis- 
tribution. The effect of this Coulomb distortion may be 
simulated in the Helm model by replacing the momen- 
tum transfer q by an effective value 15011 . In our study of 
electron scattering in isotopic chains [28j, we found that 
this correction was important because we were compar- 
ing electron-nucleus scattering observables of nuclei with 
the same proton number, and therefore they exhibited 
rather small differences. However, this correction can be 
neglected for the purposes of the present work because it 
is almost totally masked by the much larger effect from 
the changing proton number along the isotonic chain. 



III. RESULTS: N = 82, N = 50, AND N = 14 
ISOTONIC CHAINS 

We start with the discussion of the results for the N = 
82 isotonic chain where the different aspects of our study 
are described in detail. After that, we extend our study 
to the TV = 50 and N = 14 isotonic chains. 
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FIG. 1: Energy of the proton single-particle levels for icpZr, 
58°Ce, 64 6 Gd and 72 4 Hf as computed with the relativistic nu- 
clear mean field interaction G2. 
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FIG. 2: Charge densities for Jg 2 Zr, U°Ce, Gd and 7 | 4 Hf as 
a function of the radial distance to the center of the nucleus 
according to the covariant model G2 (solid lines) and to the 
fitted Helm distributions (dashed lines). 



A. N — 82 isotonic chain 

We first analyze the charge densities along the N = 82 
chain. The ordering and the energy of the different pro- 
ton single-particle levels, mainly the levels closest to the 
Fermi level, are quite important for the present study. 
This is because the corresponding single-particle wave 
functions determine, to a large extent, the shape of the 
charge density at the surface region as well as the electric 
charge form factor in the low-momentum transfer region. 
Fig. [1] displays the energy of the proton single-particle 
levels of some selected nuclei of the N = 82 isotopic 
chain. They are representative of proton deficient nuclei 
(4Q 2 Zr), stable nuclei (^Ce), proton rich-nuclei (g| 6 Gd) 
and proton drip- line nuclei ( 7 | 4 Hf). 

The more relevant proton single-particle levels in our 
analysis of the N — 82 isotonic chain are, on the one 
hand, the I39/2, I57/2 an d 2d 5 / 2 levels (which appear 
clearly separated in energy) and, on the other hand, the 
nearly degenerate l/i-ii/2) 2^3/2 and 3si/ 2 levels (which 
have a very close energy). The I/111/2 level shows energy 
gaps of about 2 and 4 MeV with respect to the 2d 5 / 2 and 
l<77/2 levels, respectively, and a gap of about 9 MeV with 
respect to the deeper 1 g$ / 2 level. This large energy gap is 
due to the magicity of the proton number Z = 50. With 
increasing mass number these relevant levels are shifted 
up in energy, roughly as a whole, retaining the same or- 
dering and approximately the same energy gaps. As a 
consequence of this level scheme, in going from the nu- 
cleus 4g 2 Zr to 58° Ce the charge densities differ basically 
by the effects of filling up the I39/2 and 1 177/2 shells, and 
in going from sf^Ce to g 46 Gd the charge densities differ 
by the occupancy the 2^5/2 shell. In these proton-rich 
isotones with mass number above A — 140, the pairing 
correlations play a non-negligible role and therefore the 
charge densities also get contributions from the I/2.11/2, 



2d 3 / 2 , and 3s!/ 2 orbitals. Finally, in the case of the drip- 
line nucleus 7 2 4 Hf, all of the mentioned single-particle 
wave functions contribute significantly to the charge den- 
sity. The differences in the charge distribution due to 
single-particle effects become evident in Fig. [2] where the 
charge densities of l 22 Zr,if°Ce, JfGd, and 7 | 4 Hf com- 
puted with the relativistic mean field model G2 are dis- 
played as functions of the radial distance. 

The equivalent Helm charge densities of these isotones, 
with parameters determined as explained in Section fll A( 
are depicted in Fig. [5] by dashed lines. As in the case of 
isotopes studied in Ref. [28j], the quantal oscillations of 
the mean-field charge densities are nicely averaged by the 
bulk part of the Helm model densities. In spite of the 
fact that the surface fall-off of the Helm densities is of 
Gaussian type, the agreement at the surface between the 
mean-field and the equivalent Helm charge distributions 
is in general satisfactory. 

The fitted parameters of the equivalent Helm densi- 
ties along the N = 82 isotopic chain, namely Rq and a 2 , 
are displayed in the two panels of Fig. [3] as functions of 
A 1 / 3 and A, respectively. The radius i?o roughly follows 
a linear trend with A 1 ^ as it can be expected from the 
increase of the total number of nucleons. The param- 
eter (T 2 , which determines the surface thickness of the 
charge density, shows a non-uniform variation with the 
mass number A, reflecting the underlying shell structure 
of the nuclei of this chain. 

In our study of isotopic chains [281 ], we found a simi- 
lar behavior of the a 2 parameter but with two important 
differences. First, the range of variation of a 2 in iso- 
topic chains is much smaller than the one exhibited by 
the N = 82 isotonic chain. Second, in the case of the Sn 
isotopes (see Fig. 6 of Ref. 28]) a 2 displays local min- 
ima for 132 Sn and 176 Sn, pointing out the magicity of the 
N = 82 and N = 126 neutron numbers which makes the 
charge densities of these isotopes more compact. In con- 
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FIG. 3: Upper panel: Mass-number dependence of the Helm 
parameter Ro predicted by the covariant mean field model 
G2 in the N = 82 isotonic chain. Lower panel: Mass-number 
dependence of the Helm parameter a 2 . The average value is 
depicted by a horizontal dashed line. 



trast, in the N = 82 isotonic chain, the kinks shown by a 2 
are rather related with the filling of the different proton 
single-particle orbitals belonging to the major shell be- 
tween Z = 50 and Z = 82. In particular, when the lgg/2 
and l<?7/2 shells are being filled, i.e., between \q*Zt and 
l^Ce, a 1 decreases almost linearly. The local minimum 
of a 2 for 5|°Ce points to some magic character of this 
nucleus. The fact that the a parameter takes the smaller 
values in the region around 53° Ce, indicates that the sur- 
face of the equivalent charge density is more abrupt at 
and around this nucleus. When the 2d§/2 level starts to 
be appreciably occupied, a 2 increases again nearly lin- 
early till g4 6 Gd, where a new kink appears. From g| 6 Gd 
to the proton drip line (^Hf), the value of a 2 continues 
to increase, but now with a smaller slope as a conse- 
quence of the higher occupancy of the l/in/2, 2^3/2, and 
3si/2 levels. Therefore, the Helm model parameter a 2 is 
extremely sensitive to the evolution of the proton shell 
structure along an isotonic chain. 

We next inspect the main properties of the differen- 
tial cross sections and electric charge form factors of the 
TV = 82 isotones. In the upper panel of Fig. @] we display 
for three representative nuclei of the N = 82 chain the 
DCS as a function of the scattering angle 8. The electron 
beam energy is 500 MeV. The DCS is computed in the 
DWBA using both the self-consistent mean-field charge 
densities obtained with the G2 model (solid lines) and 
the equivalent Helm charge densities (dashed lines). The 
square modulus of the DWBA electric charge form fac- 
tor |-Fdwba(<z)| 2 as a function of the momentum transfer 
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FIG. 4: DCS for elastic electron-nucleus scattering (upper 
panel) and square charge form factor (lower panel) in io 2 Zr, 
58°Ce, and £| 4 Hf at 500 MeV computed in DWBA. The results 
are shown both for the self-consistent mean-field densities of 
G2 (solid lines) and for the equivalent Helm distributions fit- 
ted to the G2 densities (dashed lines). In the lower panel, we 
also show by empty symbols the results obtained at 250 MeV 
using the self-consistent G2 densities. 



q = 2E sin(9/2) is shown in the lower panel of Fig. [4] 
The empty symbols in the lower panel of this figure cor- 
respond to |i 7 bwBA('?)| 2 computed at an electron beam 
energy of 250 MeV. The comparison of the results for 
|i ? DWBA(g)| 2 at E = 500 MeV and E = 250 MeV shows 
that the electric charge form factor defined in Eq. ([I} is 
largely independent of the energy of the beam in the low- 
momentum transfer domain. Therefore, the analysis of 
|-Fdwba(<7)| 2 contains the essential trends of the elastic 
electron-nucleus scattering in this regime. 

The dashed lines in the two panels of Fig. @] corre- 
spond to the DWBA result but using the equivalent Helm 
charge distributions, fitted as explained previously, in- 
stead of the self-consistent mean-field densities. One 
notes an excellent agreement at low-momentum trans- 
fers up to 1.5-2 fm -1 between the results from the orig- 
inal mean-field densities and from the equivalent Helm 
charge densities. This fact reassures one of the ability of 
the parametrized Helm distributions to describe global 
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FIG. 5: Evolution of the square modulus of the DWBA elec- 
tric charge form factor with the momentum transfer q along 
the N — 82 isotonic chain as predicted by G2 at an elec- 
tron beam energy of 500 MeV. The momentum transfer corre- 
sponding to the first inflection point for each isotone is shown 
by circles. 



trends of elastic electron-nucleus scattering at low q, as 
it was also found in Ref. [28[ for isotopes. 

In medium and heavy mass nuclei, the first oscillations 
of the DCS and of the square charge form factor com- 
puted within the DWBA usually do not show clean local 
minima but they rather show inflection points. As we can 
see in Fig. [4] this is the situation for the first oscillation 
of the DCS and of |Fdwba(q)| 2 in the N — 82 isotonic 
chain. In the absence of an explicit minimum, the first 
inflection point (IP) is the best candidate to characterize 
the relevant properties of the electric charge form factor 
at low q as we discussed in Ref. [28]. In Fig. [5] we plot 
|-Fdwba(<z)| 2 for the TV = 82 isotones in a magnified view 
around the first IP. The value of |-Fdwba(<z)| 2 at the first 
IP is depicted by circles for each nucleus. In agreement 
with earlier literature (29j , the momentum transfer at the 
first inflection point (qip) shows an inward shifting and 
the value of |-Fdwba(5ip)| 2 shows an upward trend with 
increasing mass number along the isotonic chain. 

Let us discuss possible correlations of the DWBA 
charge form factor at low-momentum transfer with the 
parameters Rq and a of the equivalent Helm charge den- 
sity, as we did in our previous analysis of isotopic chains 
[28j . If we first look at the analytical expression of the 
charge form factor predicted by the Helm model, cf. Eq. 
(HI), it suggests to use qRo and <r 2 g 2 as the natural vari- 
ables to investigate the variation of this quantity. In the 
q — > limit, Eq. ((4]) can be written as 

^(9 0) = 1 - ^q 2 (5a 2 + i? 2 ) + 0[q% (7) 

This result points towards a linear correlation with the 
mean square radius (r 2 ) of the Helm distribution (Hoj 




U , I , I , I . ^1/21 

0.3 0.4 0.5 0.6 0.7 



2 2 

° 9* 

FIG. 6: Square modulus of the electric charge form factor in 
DWBA at the first inflection point (gip) as a function of a 2 qi P 
predicted by the RMF model G2 for the N = 82 isotones. 



due to the fact that 

(rl) = l(5a*+Rt). (8) 

However, the correlation suggested by this approxima- 
tion is not fulfilled by the DWBA calculations in the 
relevant region of momentum transfers for our study. 
For this reason, we have further investigated the relation 
of |-Fdwba(<Zip)| 2 with qf P R% and er 2 g 2 P separately. We 
show in Fig. [5] the behavior of |-Fdwba(9ip)| 2 as a func- 
tion of the value of <J 2 q 2 P since it will be very instructive 
to understand the influence of the proton shell structure 
on elastic electron scattering in the isotonic chains. 

The non-uniform variations seen in Fig. [5] along the 
horizontal axis are basically due to the Helm parame- 
ter a 2 rather than to q 2 P . This is because of the fact 
that if we compare the relative change along the isotopic 
chain found in the quantities a and q\p (cf. Figs. |3] and 
[5J respectively), it is much larger in the case of the a 
parameter. Hence, the information along the horizontal 
axis of Fig. [5] is sensitive to the filling order of the single- 
particle levels contributing to the charge density at the 
surface region. 

The non-uniform variation shown by |-Fdwba(9ip)| 2 
along the vertical axis in Fig. [5] can be qualitatively un- 
derstood in terms of the single-particle contributions to 
the PWBA electric charge form factor. To this end we 
plot in Fig. [7] the contribution to the PWBA form factor 
from the individual proton orbitals: 

f nlj (q) = J dr\yj nlj (f) | 2 e** (9) 

where ip n ij (r) is the wave function of a proton level with 
quantum numbers n, I, and j. Note that @ does not 
include the occupation probability factors (i>«y) and de- 
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FIG. 7: Shell contribution to the form factor of the last oc- 
cupied levels in PWBA [see Eq. ([9])] as a function of the mo- 
mentum transfer, calculated in the nucleus 72 4 Hf with the G2 
mean field model. The shaded region indicates the range of 
observed gip values in the calculations of the form factor in 
DWBA for the N = 82 isotonic chain. 
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FIG. 8: Square modulus of the electric charge form factor in 
DWBA at the first inflection point (qip) as a function of Rq 
predicted by the RMF model G2 for the N = 82 isotones. 



generacies (2j + 1), i.e., 

FpwBA(q) = ^ ^2(2j + l)v nlj f nlj (q). (10) 

nlj 

In Fig. [7] we depict f n ij (q) for the orbitals close to the 
Fermi level in the nucleus 7 2 4 Hf, which can be consid- 
ered as representative of the level scheme of the N = 82 
isotopic chain. First, one notes that the single-particle 
contributions f n ij(q) to the PWBA electric charge form 
factor do not have the same sign in the range of mo- 
mentum transfers of interest in our analysis. Therefore, 
strong interference effects may occur among these single- 
particle contributions. In particular, we can see in Fig. [7] 
that in the region around qip the contributions from the 
lpg/2, l<?7/2i an d orbitals are negative, while the 

contributions from the 3sj/2, 2d 5 / 2 , and 2d 3 / 2 orbitals are 
positive. The PWBA form factor corresponding to the 
underlying Z = 40 core is negative. Therefore, when the 
I59/2 and l<77/ 2 orbitals are occupied — in passing from 
4Q 2 Zr to 5|°Ce — the square modulus of the PWBA form 
factor increases. When on top of this configuration, the 
2d 5 /2 orbital is filled in g| 6 Gd, the square modulus of the 
PWBA charge form factor decreases due to the positive 
sign of the contribution of this level around grp. This 
simple pattern in the uniform filling picture is slighlty 
modified due to the pairing correlations that introduce 
additional mixing with the contributions from the 2d 3 / 2 , 
3si/ 2 , and lft-n/2 orbitals. In spite of this, the simple 
PWBA description is quite useful to help us interpret 
the changes of |F DW BA(gip)| 2 from }§ 2 Zr to If Gd. The 
subsequent increase shown by |-Fdwba(<Zip)| 2 from g| 6 Gd 
to y| 4 Hf can also be understood in this schematic picture 
since the PWBA charge form factor of g| 6 Gd is globally 
negative and the additional contribution of the l/in/ 2 



orbital also is negative for q values near qip. 

The discussed theoretical results pinpoint the impor- 
tance of the filling order of the proton single-particle lev- 
els in elastic electron scattering off exotic nuclei. There- 
fore, future experiments such as those planned in the 
upgrades of the GSI and RIKEN facilities may become 
excellent probes of the shell structure of exotic nuclei and 
may also confirm or refute the theoretical predictions of 
new magic numbers. 

Regarding the relation of |Fdwba(9ip)| 2 with (jfpi? 2 ,) 
we have not found a simple behavior. In spite of this, 
we have observed that |Fdwba(9ip)| 2 and the square of 
the Helm radius Rq show a rather similar behavior as a 
function of the mass number in the isotonic chain. This 
suggests plotting |Fdwba(9ip)| 2 against Rq, which we 
do in Fig. [51 One observes a good linear correlation 
between both quantities. This correlation indicates that 
the parameter of the Helm model which measures the 
size of the bulk part of the density profile of each isotone 
governs the magnitude of the electric charge form factor 
at low momentum transfer. 

To conclude this section, we would like to note that 
some of the details of the predicted single-particle ener- 
gies, energy gaps, and filling order of the orbitals change 
to some extent if in our calculations we use other RMF 
models or Skyrme forces instead of the G2 interaction. 
In particular, this is due to the fact that we are explor- 
ing regions of the nuclear chart beyond the region where 
the parameters of these effective nuclear interactions have 
been calibrated. However, the basic conclusion to be em- 
phasized, i.e., the manifest sensitivity of some electron 
scattering observables to the proton shell structure of the 
isotones, is a robust feature that comes out regardless of 
the effective nuclear model. 



8 



70 Ca 

Z=50 



Z=28 



' 2 I',r 
if-"- 



Z=20 



Z=8 



84 Se 

Z=50 



Z=28 



Z=20 



Vm 
- w sa 



Z=8 



90 



Zr 





l J5/2 








Z=28 








Z=20 




w M 



Z=8 



100 



Sn 



2 Pl/2 



Z=28 



Z=20 



Z=8 



- 'Pin 



FIG. 9: Energy of the proton single-particle levels for 2oCa, 
IfSe, IcjZr, and so°Sn as computed with the G2 parameter set. 
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B. TV = 50 isotonic chain 

The more relevant proton single-particle orbitals for 
our study of the N = 50 chain are the I/7/2, I/5/2, 
2p 3 / 2 , 2pi/2, and lgg/2 orbitals. They cover two major 
shells between Ca and Sn. This set of proton energy lev- 
els computed with the G2 parametrization is displayed 
in Fig. [9] for some selected isotones. We can see that 
these levels move up in energy, roughly as a whole, when 
the mass number increases in going from proton-deficient 
nuclei ( 2 (jCa) to stable nuclei (ffSe, foZr) and to proton 
drip-line nuclei (j^Sn). 

The parameters i?o an d cr 2 of the Helm model distri- 
butions fitted to the mean-field charge densities of the 
N = 50 isotones are displayed in Fig. [TO] The global 
features are similar to the case of the N = 82 chain. The 
Ro parameter, which represents the effective location of 
the surface of the nucleus, approximately follows a lin- 
ear trend with A 1 / 3 . The mass- number dependence of 
<j again displays a non-uniform trend, originated by the 
filling of the different proton single-particle orbitals. We 
see that a 2 decreases in filling the I/7/2 shell from 2oCa 
to 2 |Ni, it remains roughly constant when the I/5/2 level 
is being filled up to ||Se, it increases when the 2y> 3 / 2 
and 2p\/2 shells are occupied till foZr, and it then de- 
creases until the proton drip-line nucleus 5o°Sn is reached 
by filling the l<?g/ 2 level. Therefore, the more abrupt 
(smaller a) equivalent charge densities predicted by the 
G2 model in the N = 50 isotonic chain correspond to 
nuclei between the doubly-magic, proton-deficient jfNi 
nucleus and the more stable 34 Se nucleus, where mainly 
the I/5/2 shell has been filled. In these nuclei, the oc- 
cupancy of the 2p 3 / 2 , 2p!/ 2 , and I59/2 levels due to the 
pairing correlations is rather small. 

The square modulus of the DWB A electric charge form 
factor for an electron beam energy of 500 MeV is dis- 
played against <J 2 q\ v in Fig. [TT] As in the case of the 



FIG. 10: Upper panel: Mass-number dependence of the Helm 
parameter 7?o predicted by the covariant mean field model 
G2 in the N = 50 isotonic chain. Lower panel: Mass-number 
dependence of the Helm parameter a 2 . The average value is 
depicted by a horizontal dashed line. 



N = 82 isotones, the behavior of <J 2 q 2 p is dominated by 
the Helm parameter a. This is because the relative vari- 
ation of a 2 (see Fig. H0| is much larger than the relative 
variation of along the isotonic chain. The change 
of |Fdwba(<Zip)| 2 along the N — 50 chain shows, glob- 
ally, an increasing trend with the mass number. We 
can appreciate in Fig. [TTJ that although the variation 
of I-Fdwba^ip)! 2 is almost linear when a specific pro- 




FIG. 11: Square modulus of the electric charge form factor in 
DWBA at the first inflection point (gip) as a function of a 2 q 2 P 
predicted by the RMF model G2 for the N = 50 isotones. 
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FIG. 12: Square modulus of the electric charge form factor in 
DWBA at the first inflection point (gip) as a function of Rq 
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ton orbital is being occupied, drastic changes of slope 
take place when a new shell starts to be significantly 
occupied. Recalling the simplified PWBA picture, cf. 
Eqs. © and (|10[) . we find that in the region of q values 
around qip the contribution to the electric charge form 
factor from the 1/ and lg orbitals is negative, while the 
contribution from the 2p orbitals is positive. This fact 
is consistent with the behavior shown by |-FbwBA(<Zip)| 2 
in Fig. [11] That is, |-Fdwba(<Zip)| 2 increases in pass- 
ing from ™Ca to ||Se, basically due to the filling of the FIG. 13: Energy of the neutron (upper panel) and proton 
I/7/2 an d I/5/2 shells, and then its value is practically 
quenched up to 4(JZr because the 2p 3 / 2 and 2p±/2 orbitals 
contribute with opposite sign to the If orbitals. When 
the lgg/2 level is appreciably occupied in approaching 
the proton drip line, the value of -Fdwba^ip)! 2 increases 
again with a nearly constant rate. Finally, in Fig. [T2l we 
see that -Fdwba^ip)! 2 of the N = 50 isotones shows 
a good linear correlation with the square of the Helm 
parameter Rq. 
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(lower panel) single-particle levels for g O, ioNe, i2>Mg, 14S1, 
fgS, fgAr and loCa as computed with the G2 parameter set. 



C. TV = 14 isotonic chain 

In this section we discuss the lightest isotonic chain an- 
alyzed in our work. Although the present findings are to 
be taken with some reservations because the mean-field 
approach is not best suited for light-mass exotic nuclei, 
we note that similar general trends to those observed 
in the heavier-mass isotonic chains also appear in the 
N = 14 chain. 

In Fig. Q2] we display the neutron (upper panel) and 
proton (lower panel) single-particle levels computed with 
the G2 interaction for the N = 14 isotonic chain, from the 
very proton-deficient nucleus | 2 to the very proton-rich 
nucleus j^Ca. As expected, the neutron single-particle 



levels become more bound with increasing mass number. 
In addition to the prominent energy gap at N — 8 seen in 
the whole chain, it may be noticed that the ld 5 / 2 neutron 
level becomes progressively more isolated when the mass 
number increases, which points to some magic character 
of the neutron number N = 14 in the calculation with the 
G2 model. This magic character is confirmed by the van- 
ishing neutron pairing gap found in our calculation from 
fjMg to 2oCa. It may be observed that the G2 model 
also predicts a slightly magic character of the neutron 
number N = 16 towards the neutron drip line. Actu- 
ally, we see in Fig. Q2] that the relatively magic trend of 
N = 14 increases from the proton-poor side (| 2 0) to the 
proton-rich side (toCa) of the chain, while the somewhat 
magic trend of N = 16 decreases from | 2 to 20 ^a. 

The more relevant proton single-particle orbitals for 
our study of the N = 14 chain belong to the s-d major 
shell. The energy levels of this proton major shell (see 
the lower panel of Fig. I13[) lie approximately at the same 
energy for all the nuclei from 22 O to 20 Ca, with roughly 
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FIG. 14: Evolution of the square modulus of the DWBA elec- 
tric charge form factor with the momentum transfer q along 
the N = 14 isotonic chain as predicted by G2 at an elec- 
tron beam energy of 500 MeV. The momentum transfer cor- 
responding to the first minimum for each isotone is shown by 
circles. 



~0 0.5 1 1.5 2 2.5 3 

q (fm ) 

FIG. 15: Shell contribution to the form factor of the last 
occupied levels in PWBA [see Eq. |(9|] as a function of the 
momentum transfer. The shaded region indicates the range 
of observed q m i n values in the calculations of the form factor 
in DWBA for the N = 14 isotonic chain. 




constant energy gaps. The ld§/2 and 2si/ 2 proton lev- 
els exhibit a considerable energy gap between them in 
this isotonic chain according to the predictions of the G2 
model. It is to be mentioned that due to the pairing cor- 
relations, the proton levels I/7/2 and I/5/2 (the latter is 
not displayed in Fig. IT3|) also play some role in our cal- 
culation of the mean-field charge densities. These levels 
simulate to a certain extent the effect of the continuum 
due to their quasi-bound character owing to the Coulomb 
and centrifugal barriers [5J]. 

In Fig. [14] we display a magnified view of \FuwBA(q)\ 2 
against the momentum transfer for the TV = 14 isotonic 
chain. We see that in this chain of lower mass, the in- 
flection point that was found after the first oscillation 
of the charge form factor in the heavier chains N = 50 
and N = 82 becomes a clearly well defined local min- 
imum. Thus, for the discussion of the N = 14 chain 
we focus on the properties of |-Fdwba(<7)| 2 at its first 
minimum (q m i n ). The value of |i 7 bwBA(9min)| 2 (shown 
by the circles in Fig. [14]) increases when the value of 
the momentum transfer at the first minimum decreases, 
i.e., |i 7 bwBA(9min)| 2 grows with increasing mass number 
in the chain. Though the increase of |-FbwBA(9mm)| 2 m 
Fig-HUis roughly linear with q m in, one notes a kink at the 
point corresponding to the 2 |Si nucleus. As we can realize 
from Fig. [12] in the schematic PWBA picture, this kink 
is originated by cancellation effects between the oppo- 
site contributions to the charge form factor around q m i n 
coming from the single-particle wave function of the Id 
proton level (negative contribution) and of the 2s proton 
level (positive contribution). 

The upper panel of Fig. [!()] shows that the parameter 
Rq of the equivalent Helm charge densities displays, as 
in the heavier isotonic chains, an overall linear increas- 



ing trend with A 1 / 3 . In turn, the variation of the Helm 
parameter a reflects the underlying shell structure of the 
mean-field charge densities. In the lower panel of Fig. 1161 
we see that the value of a 1 remains almost constant be- 
tween g 2 and f|Si when mainly the lc^/2 shell is being 
filled. From 3 gS on, the 2si/ 2 and ld 3 / 2 levels start to 
be appreciably occupied and a 2 starts increasing almost 
linearly with A till the proton-drip line nucleus 2oCa. 
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The influence of the discussed proton shell structure 
on the electric charge form factor at the first minimum 
for the N = 14 isotones is obvious in Fig. [TT1 which 
displays \F DWBA (q min )\ 2 against the value of cr 2 g 2 lin . Fi- 
nally, if we analyze the variation of -FbwBA(9min)| 2 with 
the Helm parameter Rq, a correlation is found between 
both quantities (cf. Fig. [T8]) . though now this correlation 
is less linear than in the heavier isotonic chains TV — 82 
and N = 50 (cf. Figs. M and 
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FIG. 18: Square modulus of the electric charge form factor 
in DWBA at the first minimum (g m in) as a function of Rq 
predicted by the RMF model G2 for the N — 14 isotones. 



IV. SUMMARY AND CONCLUSIONS 

In this work we have explored some of the information 
about nuclear structure that can be obtained from elas- 
tic electron scattering in isotones. Due to the fact that 
the number of protons changes along an isotonic chain, 
this study primarily probes the effect of the different pro- 
ton single- particle shells on the elastic electron scattering 
observables. 

We have computed the DWBA differential cross sec- 
tion at an electron beam energy of 500 MeV using the 
charge densities calculated self-consistently with the co- 
variant nuclear mean- field model G2 [HEl. The electric 
charge form factor has been obtained by taking the ra- 
tio of the DWBA differential cross section with the exact 
(DWBA) Mott differential cross section. The so-defined 
electric charge form factor is practically independent of 
the electron beam energy in the low-momentum transfer 
regime. 

We have paid special attention to the electric charge 
form factor taken at the momentum transfer qip of the 
first inflection point, for the N = 50 and N — 82 chains, 
or at the momentum transfer g m i n of the first minimum 
for the N = 14 chain. In agreement with earlier litera- 
ture, we have found that the values of gip ((? m i n ) shift 
inwards (i.e., they become smaller) and that the val- 
ues of |F DW ba(<7ip)| 2 (l-FbwBA(gmin)l 2 ) increase when the 
atomic number increases from the neutron drip line to the 
proton drip line of the isotonic chain. 

The results reveal that along an isotonic chain the 
DWBA electric charge form factor of each nucleus is ex- 
tremely sensitive to the underlying proton shell struc- 
ture. In the simpler PWBA picture, we have found that 
a particular proton shell may contribute to the electric 
charge form factor with positive or negative sign in the 
momentum transfer region of interest in our analysis. 
The contributions from levels with more than one radial 
node and small orbital angular momenta have opposite 
sign to the contributions from levels without radial nodes 
and large angular momenta. As a consequence, cancella- 
tion effects can appear when the different proton shells 
are successively occupied. Although this description is 
to some extent masked by pairing correlations, a similar 
situation is found in the DWBA calculations of the elec- 
tric charge form factor. Therefore, the rate of change of 
|-Fdwba(9ip)| 2 or |F DWB A(gmin)| 2 with the proton num- 
ber along an isotonic chain may vary substantially when 
a new single-particle orbital enters the nucleus. This sug- 
gests that scattering experiments performed on isotones 
can be effective probes of the proton nuclear shell struc- 
ture of stable and unstable nuclides. 

To investigate the dependence of the first inflection 
point or first minimum of the electric charge form factor 
with the basic properties of the charge distribution, we 
have parametrized the moan-field densities by the Helm 
model. We have found that the Helm parameter Ro, 
which measures the mean position of the surface of the 
charge density, increases with the mass number following 
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roughly an A 1 / 3 law. The Helm parameter a, which is 
related with the surface thickness of the charge distri- 
bution, encodes information from the underlying proton 
shell structure and, as a consequence, it does not show 
a definite trend with the mass number. We have also 
found that at low-momentum tranfers the square modu- 
lus of the DWBA electric charge form factor is accurately 
reproduced if the mean-field densities are replaced by the 
fitted Helm densities. 

In a previous work [28| , we noted that correlations be- 
tween the DWBA electric charge form factor at the first 
inflection point (or minimum) and the parameters of the 
Helm densities can provide information about global fea- 
tures of electron scattering along isotopic chains. In the 
isotonic chains analyzed in the present work, we find that 
|-Fdwba(<7ip)| 2 (l^bwBA^min)! 2 ) shows a rather good lin- 
ear correlation with the Helm parameter R^, specially 
in the heavier isotonic chains, while there is no regular 
behavior with the Helm parameter <r 2 because of its de- 
pendence on the last occupied proton orbitals. It should 
be pointed out that shell effects encoded in the DWBA 
electric charge form factor are magnified if |FbwBA(<Zip)| 2 
(|-FDWBA(<Zmin)| 2 ) is plotted against <j 2 qf P (cr 2 ? 2 ^) along 



an isotonic chain, providing interesting insights into the 
proton shell structure of the nuclei of the chain. 

In summary, we can conclude that electron scattering 
on unstable nuclei that is intended to be performed in 
next-generation facilities [l5|-[22j , may become a very use- 
ful tool to obtain relevant information about the single- 
particle shell structure of exotic nuclei. 
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